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Abstract 


Hus document contains graphical data ncccssai) lot the pichminai) design of ballistic 
missions to Jupiter Con touts of launch energx lequirements. as well as man> othei 
I.iuncli and Jupilci amval paiamctcis. ate piesenied in launch date airival date space 
toi j!I launch oppoitunities fioni 1**85 tluough 2005. In addition, an extensive text is 
included which explains mission design methods, tiom launch window development to 
Jupitei piobe and oihitei amval design, utilising the giaplucal data m this volume as 
well as mimeious equations lelating vauous paiamctcis. Fins is the fust of a planned 
seiies M mission design documents which w ill appl> to all planets and some othei bodies 
in the solai s\stem. 



Preface 


This publication is one of a series of volumes devoted to interplanetary trajectories of 
different types. Volume 1, of which the present publication is Part 3, describes ballistic 
trajectories. Parts 1.2, and 4, which will be published in 1983. will treat ballistic trajec¬ 
tories to Venus. Mars, and Saturn, respectively. 
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I. Introduction 

The purpose of tins series of Mission Design Handbooks is 
to provide trajectory designers and mission planners with 
graphical trajectory information, sufficient for preliminary 
interplanetary mission design and evaluation. In most respects 
the series is a continuation of the previous three volumes of 
the Mission Design Data, TM 33-73 ft (Ref. 1) and its predeces¬ 
sors (e.g., Ref. 2); it extends their coverage to departures 
through the year 3.003 A.D. 


but certainly not unrealistic using future orbital assembly 
tt'rhniqiifs. 

A separate series of volumes (Ret. 3) is being published 
concurrently to provide purely geometrical (i.e., trajectory- 
independent) data oi planetary positions and vfcwing/nncnta- 
tion angles, experienced by a spacecraft in the vicinity of these 
planetary bodies. The data cover the time span through 2020 
A.D.. in order to allow sufficient mission duration time lor all 
Earth departures, up to 2005 A.D. 


The entire series is planned as a sequence of volumes, each 
describing a distinct mission mode as follows: 

Volume 1* Ballistic (i.e.. unpowered) transfers be¬ 
tween Earth and a planet, consisting of 
one leg trajectory arcs. For Venus and 
Mars missions the planet-to-Earth return 
trajectory data are also provided. 

Volume II: Gravity-Assist (G'A) trajectory transfers, 
comprising from two to four ballistic 
interplanetary legs, connected hy suc¬ 
cessive planetary swingbys. 

Volume III: Delta-V-EGA (AVEGA) transfer trajec¬ 
tories utilizing an impulsive deep-space 
phasing and shaping burn, followed hy a 
return to Earth for a G/A swingbv maneu¬ 
ver taking the spacecraft (SO to the 
eventual target planet. 

Each volume consists of several parts, desetibing Irujecto: 
opportunities for missions toward specific taigot oi swingbv 
bodies. 


This Volume 1. Pait 3 of the series is devoted to ballistic 
transfers between Earth and Jupiter. It ucserihes tiajectorics 
taking from 1 to 5 yeais of flight time for the 20 successive 
mission opportunities, departing Earth in the following years: 
l‘>K5. D>8b, D>8“\ 1 4 >8S. 1 4 >S 4 >, DK>(). 1 4 > 4 >I. 1 4 > 4 >2'3. 1<*I3 4. 
1004 5. l‘)%, 1*W\ loos.DK> 4 >. 2000. 2001. 2002. 2003 4. 
2004 5. and 2005 t« 


Individual variables picsented herein arc described in detail 
in subsequent sections and summarized again in Section IV’. 
Suffice it to say here that all the data aie presented in sets of 
11 contour plots each, displayed on the launch date anival 
date space foi each opportunity, kequired departure eneigy 
t\ % departure asvmptoK decimation and light ascension, 
amval 1^ and its cquatonul directions, as well > Sun and 
Eat th direction angles with i expect to the depart me auixal 
asymptotes. are presented. 


The geometric data are presented in graphical form and 
consist of 20 quantities, combined .. * > eight plots for each 
calendar year and each target planet. Tiu graphs display equa¬ 
torial declination and right ascension of Earth and Sun (plan- 
etoccntric). as well as those of the target planet (geocentric); 
heliocentric (ecliptic) longitude of the planet, its heliocentric 
and geocentric distance; cone angles of Earth and Canopus, 
clock angle of Earth (when Sun/Canopus-oricntcd); Earth-Sun- 
planel. as well as Sun-Earth-planet angles; and finally, rise and 
set times for six deep-space tracking stations assuming a h-deg 


luni/on mask. This information is similar to that in the second 
part of each of the volumes previously published (Ref. 1). 


II. Computational Algorithms 

A. General Description 

The plots >or the entire series were computer-generated. A 
minimum of editorial and graphic support was postulated from 
the outset in an effort to reduce cost. 


A number of computer programs were created and/or modi¬ 
fied to suit the needs of the Handbook production. 

The computing effort involved the genet at ion ot arrays ot 
transfer trajectory arcs connecting departure and arrival 
plane is on a large number ot suitable dates at each body. 
Algorithms (computational models) to solve this problem can 
vary greatly as to their complexity, cost ot data generated, and 
resulting data accuracy. In light of these conrideiutions, the 
choice of methods used in this etfort has been assessed. 


B. Two-Body Conic Transfer 

I* !ch depaitute annal date combination tepiesenis a unique 
tiaiiNtoi tiatectoiy between two ,itied bodies, it the uu.ubei 
ot (evolutions ot the space*.utt about the piimatv (c.c . the 
Sun) is spi\ lied I he I ambeit Iheoieni piovuies a suitable 
tiamewoik t.»i the imputation ot siiji tiaie*.tones, but 
only it leslnctcd two bodv e>iik motion ptevaiU 


It should be •• led that pait« »t the lautuf' 'pace co\ ..d 
may tequire launchet eneicies not ptesently r-'NJ) available. 


Restricted two-bodv nioiion implies that the dynamical 
system consists ot onlv two bodies, one ot which, the pmiiaiv. 
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is so much more massive than the other, that a!! ol the system’s 
gravitational attraction may he .ssumed as concentrated at a 
point the center ol that primary body. I he secondary body 
of negligible mass (c.g., the spacecraft) then moves in Kep- 
lerian (conic) orbits about the primary (e.g.. the Sun) in such 
a way that the center of the primary is located at one of the 
foci of the conic (an ellipse, paiabola, or hyperbola). 

The Lambert Theorem states that given a value of the; *uvi- 
tationa! parameter n (also known as UM) for the central body, 
the time of llight between two arbitrary points in space, 
R 1 and R 2 . is a function of only three independent variables: 
the sum of the distances of the two points front the focus, 
|R,I + IR 2 I, the distance between the two points C= |R 2 - Rj 1. 
and the semimajor axis. a. of the conic orbital flight path 
between them (Fig. I.) 



Detailed algorithm descriptions of the Lambert method, in* 
eluding nccessau branching and singularity precautions, are 
presented in numerous publications, e.g.. Refs. 2 and 4. The 
computations result in a set of conic classical elements (d. t\ 
il. cj. n,) and the transfer angle. Au l2 . or two equivalent 
spacecraft heliocentric velocity vectors. one at depar¬ 

ture. the other at the at rival planet. Suhtiaction of th»* appro 
priate planetary heliocentric velocity vector. V*. ........ at 

the two corresponding times from each of these two space¬ 
craft '.olocity vectors results in a pair of planetocentric velocity 
states “at infinity ' with respect to each planet (Fig. 2) 


v - V - V 

h S'C a II n f* /. A .V hit 


(i) 


where < = I and 2 refer to positions at departure and arrival, 
respectively. The scalar of this I* vector is also referred to as 


Fig. 1. The Lambert problem geometry 


the hyperbolic excess velocity, “^infinity" or simply “speed” 
(e.g.. Ref. 4). The represents the velocity of the spacecraft 
at a great distance from the planet (where its gravitationa 1 
attraction is practically negligible). It is attained when the 
spacecraft has climbed away from the departure planet, fol¬ 
lowing injection at velocity V.: 



or before it starts its fall into the arrival planet’s gravity well, 
where it eventually reaches a closest approach (periapse, C7.4) 
velocity L,,. 



/ IK Hi 





Fig. 2. Departure geometry and velocity vector diagram 
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1 ho variables r t and r p refer to the departure injection and 
arrival periapse planetocentric radii, respectively. Values for 
the gravitational parameter /i (or GM) are given in subsequent 
Section V on constants. 

The V w vectors, computed by the Lambert method, repre¬ 
sent a body center to body center transicr. They can, however, 
be translated parallel to themselves at either body without 
excessive error due to the offset, and a great variety cf re .istic 
departure and arrival trajectories may thus be constructed 
through their use. to be discussed later, fhe magnitude and 
direction of l’„ as well as the angles that this vector forms 
with the Sun and harth direction vectors at each terminus, are 
required for these mission design exercises. 

Misv* * tin* relatively small terrestrial planets are suited 
to 1.' analyzed by the Lambcrt method, as the problem can he 
adequately represented by the restricted two-body formula¬ 
tion, resulting in flight time errors of less than I day an 
accuracy that cannot even be read from the contour plots 
presented in this document. 

C. Pseudostate Method 

Actual precision interplanetai transfer traiectories. espe¬ 
cially those involving the giant outer planets, do noticeably 
violate the assumptions inherent in the Lambert Theorem. 
The restricted two-body problem, on which that theorem is 
based, is supposed to describe the conic motion of a massless 
secondary (i.e., the spacecraft) about the point mass of a pri¬ 
mary attractive body (i.e.. the Sun), both objects being placed 
in an otherwise empty Universe. In reality, the gravitational 
a-'ruction ot either departure or target body may significantly 
ab-.r the entire transfer trajectory. 

Numerical Vhody trajectory integration could be called 
upon to represent the true physical mode! tor the laws of 
motion, but would be too mostly, mnsidcinip the numbei ot 
complete liajes.tones requited to fully watch and desutbe a 
given misMon oppoi(unity 

The pseudostate theory, lust introduced by S. Wilson 
(Kef. 5) and modified to solve the three-body Lambert prob¬ 
lem by l). V Byrnes (Ret. t>). represents an extremely usetul 
improvement over the standard Lambert solution, Foi the 
giant planet missions, it can correc about *>5 percent ot the 
three-body errors incurred, e g., up to 30 days in flight tune on 
a Jupitci-bound tourney 


Pseudostate theory is based on the assumption that for 
modest gravitational perturbations the spacecraft conic mot.on 
about the primary and the pseudo-conic displacement due to a 
third body may be superimposed, if certain rules arc followed. 

The method, as applied to transfer trajectory generation, 
does not provide a flight path only its end states. It solves the 
original Lambert problem, however, not between the true 
planetary positions themselves, but instead, between two com¬ 
puted “pseudostates." These are obtained by iteration on two 
displacement vectors off the planetary ephemeris positions on 
the dates *»f departure and arrival. By a suitable superposition 
with a planetocentric rectilinear impact hypeibola and a 
constant-velocity, "zero gravity." sweepback at each end of 
the Lambertian conic (see Fig. 3). a satisfactory match is 
obtained. 

Of the five arcs involved in the iteration, the last three 
(towaids and at Jupiter) act over the lull flight time. A7 J2 . 
and represent: 

(1) Conic heliocentric motion between the two pseudo¬ 
states R ( and R* (capital R is used here for all helit 
centric positions). 

(2) The transformation of R* to a planetocentric position, 

r 2 (lower case r is used for planetocentric positions), 
performed in the usual manner is followed by a 
“constant-velocity" sweepback in ume to a point 
r* * r 2 - \„ 2 X correcting the planetocen 4 im¬ 

position r 2 to what it would have been at T,. had 
there been no solar attraction during ^7^, and 
finally. 

(3) The planetocentric rectilinear incoming hyperbola, 
characteu/ed by incoming L-infinity V^, a ladial 


/ \ 

S \ 

/ NO GRAVITY \ 



Fig. 3. Pseud os tata transfer gtomatry 
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nl'inpt iTnnml nnH n trin tinin \ 7" from r* 

• “■ p* • J* • ■■■< j “V »« MtlkA U ftl ft|< MHlii M 1 | 2 •••*••« * J 

to perijove. which can be satisfied by iteration on the 
r -magnitude and ;hus also on R*. 

* m 


This last aspect provides for a great simplification of the 
formulation as the R? end-point locus now moves only along 
the V, || rl vector direction. The resulting reduction in com¬ 
puting cost *s significant, and the equivalence to the Lambert 
point-to-point conic transfer model is attractive. 


The first two segments of the transfer associated with the 
departure planet may be treated in a like manner. If the planet 
is harth. the pseudostate correction may be disregarded tie . 
Rj = R f ). or else the duration ot Earth's perturbative effect 
may be reduced to a fraction of AT^. It can also be set to 
equal a fixed quantity, e g., = 20 days. The latter value 

was in fact used at the Earth's side ot the transt'c* m the 
data generation process for this document 


The rectilinear pseudostatc method thus involves an itera 
five procedure, utilizing the standard Lambert algorithm to 
obtain a starting set of values for l ^ at each end of the trans¬ 
fer arc This first guess is then improved by allowing the 
planetocentric pseudostate position vector r* to he scaled up 
and down, using a suitable partial at either bodv. such that 
AT Rh: the rime required to fall along the rectilinear hypei- 

bola through tJ (the sum of the sweepback distance l' m 
X AT and the planetocentric distance ir |). equal thegravita 

I f 

tionai perturbation duration. A7j Both t g and V along 
which the rectilinear fall occuis, are continuously leset utiliz- 
" the latest values of magnitude and direction of at each 
end. i. ot the new Lambert trariste* .uc. as the iteration pro 
grosses. The proceduie conveiges rap’.dly as the hyporholic trip 
time discrepancy. AAV - A'T Rt ^ - AT. tails below a preset 
small tolerance 


Once the V. vectors at each planet are converged upon, 
the desired output variables can he generated and contour 
plotted by existing standard algorithms. 


III. Tre.lectory Characteristics 
A. Mission Space 

All realistic launch and inaction vehicles are energy-limited 
and impose very stringent constrains on the interplan.(a*y 
mission selection process. Only those transfer opp. 
which occur near the times of a minimum Larth *•* 
energy requ.rement are thus of practical interest. H» . 
side of such an optimal date, departtite energy uicre.t *s. first 
slowly. foiKiwed by a rapid increase, thus requiring either a 
greatei launch capability, or alternatively a lower allowable 
pavload mass. A "launch period." measured hi days m even 


ii'iinlr 

■■ wva., 1 . 
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capability of a given launch/injection vehicle must equal or 
exceed the departure energy requirement for a specified 
payload weight. 


In the course of time these minimum departure energy 
opportunities do recur regularly, at “synodic period * intervals, 
reflecting a repetition of the relative angular geometry of the 
two planets. If and cj 2 are the orbital angular rates of the 
inner and outer of the two planets, respectively, moving about 
the Sun in circular orbits, then the mutual configuiation of 
the two bodies changes at the following rate: 

to = td - gj, . rad/s (4) 

I I m 

If a period of revolution. P % is defined as 



then 

J_ _ J_ 

n p p (b) 

s r i : 


where P s , the synodic period, is the period of planetary 
geometry recurrence, while P x and P 2 are the orbital "sidereal 
(i.e., inertial) periods" of the inner (faster) and the outer 
(slower) planet considered, respectively. 

Since planetary orbits are neither exactly circular nor cu- 
planar. launch opportur :, ies do not repeat exactly, some 
years being better than others in energy requirements or in 
other parameters. A complete repeat of trajectory character¬ 
istics occurs only when exactly the same orbital geometry of 
departure and arrival body recurs. For negligibly perturbed 
planets approximately identical inertial positions in space at 
departure and arrival imply near-recurren.ce of transfer tra- 
jectorv characteristics. Such events cjn rigorously he assessed 
only for nearly resonant nonprecessing planet;**) orbits, i.e., 
tor those whose periods can he related in terms of integer 
fractions. For instance, if five revolutions of one bodv cor- 
respond to three revolutions of the i ther, that time interval 
would constitute the “period of repeated characteristics." 
Near-integer ratios provide nearly petitive configurations 
with respect to the lines ot apsides and nodes. The Earth- 
relative synodic perunl ot Jupiter is >:84 days, i.e., about 
13 months. LjcIi cycle ot 11 consecutive Jovian mission 
opportunities amounts to 4.'X 1 7 d/ys and is nearly repetitive, 
driven by Jupiter’s sidereal period of 11 SrOO veals »r 4332.b 
days It is obvious that lor an identical mutual angular geom 
etiv Jupiter would be tumid beyond its nurtial positrm in the 
previous vV«,ie hv I days worth ot motion (4 vs dee I wni!«* 
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Faith would have completed 4JK7.7Alh.v25 - 12.01 J levolu 
Hons, being ahead ol the old maik by the same angulut aiiioiiut 
as Jupitei 

A variety of consideiations fore** the realist! launch period 
not to occui at the minimum energy combi'* non ol depai- 
turc and arrival dates, launch vehicle readiness status, proce¬ 
dure slippage, wearier anomalies, multiple launch strategies, 
arrival characteristics all cause the launch or. more generally, 
the departure period to be extended ovei a number of days or 
weeks and not necessarily centered on the minimum energy 
date. 


For this document, a 120-day departure date coverage span 
was selected, primarily in order to encompass launch energy 
requirements of up to a ~ 200 km*/s 2 contour, where 
C\ = rj i.c., twice the injection energy per unit mass. 


/*./ r 1-/2. Arrival date coverage was set at 1600 days to dis¬ 
play missions tiom 1- to 5-years flight time. 

The matrix of departure and arrival dates to he presented 
comprises the “mission space” for each departure opportunity. 


B. Transfer Trajectory 


As previously stated, each pair of departure/arrival dates 
specifies a unique transfer trajectory. Fach such point in the 
mission space liar* associated with it an array ot descriptive 
variables. Departure cnctgy. characterized by is by tar the 
most significant among these parameters. It increases towards 
the '*dges of the mission space, hut if also experiences a dru- 
nulic rise along a “ridge ” passing diagonally from lower left 
to upper right across the mission space (Fig 4). This distur- 
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banco is associated with all diametric, i.e., near-180-deg, 
transfer trajectories (Fig. 5). 

In 3-uimensional space the fact that all planetary orbits are 
not strictly coplanar causes such diametric transfer arcs to 
require high ecliptic inclinations, culminating in a polar flight 
path for an exact 180-deg ecliptic longitude increment be¬ 
tween departure and arrival points. The reason for this 
behavior is, as shown in Fig. 5, that the Sun and both trajec¬ 
tory end points must lie in a single plane, while they are also 
lining up along the same diameter across the ecliptic. The 
slightest target planet orbital inclination causes a deviation 

= TRUE TRANSFER 
ANGLE 

= TRANSFER ANGLE 
POLAR IN ECLIPTIC PLANE 



Fig. 5. Effect of transfer angle upon inclination of 

trajectory arc 



Fig. 6. Nodal transfer geometry 


out of the ecliptic and forces a polar 180-deg transfer, in order 
to pick up the target’s vertical out-of-plane displacement. 

The obvious sole exception to this rule is the nodal transfer 
mission, where depaiture occurs at one node of the target 
planet orbit plane with the ecliptic, whereas arrival occurs at 
the opposite such node. In these special cases, which recur 
every half of the repeatability cycle, discussed in the preceding 
paragraph, the transfer trajectory plane is indeterminate and 
may as well lie in the departure planet’s orbit plane, thus 
requiring a hsser departure energy (Fig. 6). The opposite strat¬ 
egy (i.e., a transfer in the arrival planet's orbit plane) may be 
preferred if arrival energy, V.* 2 . is to be minimized (Fig. 7). 

It should be noted that nodal transfers, being associate*, 
with a specific Jupiter arrival date, may show up in the data on 
several consecutive Jovian mission opportunity graphs (e.g., 
1991-1994/5), at points corresponding to the particularnodal 
arrival date. Their mission space position moves, from oppor¬ 
tunity to opportunity, along the 180-deg transfer ridge, by 
gradually sliding towards shorter trip times and earlier relative 
departure dates. Only one of these opportunities would occur 
at or near the minimum departure energy or the minimum 
arrival date, which require a near-perihelion to near¬ 
aphelion transfer trajectory. These pseudo-Hohmann nodal 
transfer opportunities provide significant energy advantages, 
but represent singularities, i.e.. single-time-point missions, 
with extremely high error sensitivities. Present-day mission 
planning does not allow single fixed-time departure strategies; 
however, future operations modes, e.g.. space station “on- 
time" launch, or alternately Earth gravity assist (repeated) 
encounter at a specific time, may allow the advantages of a 
nodal transfer to be utilized in full. 

The 180-deg transfer ridge subdivides the mission space 
into two basic regions: the Type I trajectory space below the 
ridge, exhibiting less th:n 180-deg transfer arcs, and the 
Type II space whose transfers are longer than 180 deg. In 
general the first type also provides shorter trip times. 

Trajectories of both types are further subdivided in two 
parts-Classes 1 and 2. These are separated, generally hori¬ 
zontally, by a boundary representing the locus of lowest 
C 3 energy for each departure date. Classes separate longer 
duration missions from shorter ones within each type. Type I, 
Class I missions could thus be preferred because of their 
shorter trip times. 

Transfer energies become extremely high for very short 
trip times, infinite if launch date equals arrival date, and of 
course, meaningless for negative trip times. 

The reason that high-inclination transfers, as found along 
the ridge, also require such high energy expenditures at depar- 
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Fig. 7. Mission space with nodal transfer 


ture is that the spacecraft velocity vector due to the Earth's 
orbital velocity must be rotated through large angles out of the 
ecliptic in addition to the need to acquire the required transfer 
trajectory energy. The value of C 3 on the ridge is large but 
finite; its saddle point minimum value occurs for a pseudo- 
Hohmann (i.c., perihelion to aphelion) polar transfer, requiring 

C, = I’? ( 7 —^ + l| * 2400 km 2 /s 2 (7) 

h \ 1+a P / 

where \\ = 2 < >.7fth km/s, the Earth's heliocentric orbital 
velocity, and a p = 5.2 Al*. Jupiter's (the arrival planet's) semi- 
majoi axis B> a similar estimate, it can he shown that for a 
true nodal pscudo-riohmann tiansfei. the minimum eneigy 
uv 1 d would reduce to 


C = u 2 

* NODAL h: 


This is the lowest value of C 3 required to fly from Earth to 
Jupiter, assuming circular planetary orbits. 

Arrival P-infinity, V«o x , is at its lowest when the transfer 
trajectory is near-coplanar and tangential to the target planet 
orhit at arrival. 

Both C 3 and near the ridge can be significantly low¬ 
ered if deep-space deterministic maneuvers are introduced into 
the mission. Tne “broken-plane" maneuvers are a category of 
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such ridge-counteracting measures, which can nearly eliminate 
all vestiges ot the near-18U-deg transfer difficulties. 

The basic principle employed in broken-plane transfers is to 
avoid high ecliptic inclinations of the trajectory by performing 
a plane change maneuver in the general vicinity of the halfway 
point, such that it would correct the spacecraft's aim toward 
the target planet’s out-of-ecliptic position (Fig. 8). 

Graphical data can be presented for this type of mission, 
but it requires an optimization of the sum of critical AV 
expenditures. The decision on which AKs should be included 
must be based on some knowledge of overall staging and arrival 
intentions, e.g., departure injection and arrival orbit insertion 
vehicle capabilities and geometric constraints or objectives 
contemplated. As an illustration, a sketch of resulting contours 
of C 3 is shown in Fig. 9 for a typical broken-plane oppor¬ 
tunity represented as a narrow strip covering the ridge area 
(Class 2 of Type I and Class 1 of Type II) on a nominal 1985 
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Fig. 8. Broken-plane transfer geometry 
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launch/arrival date C 3 contour plot. The deep-space maneuver 
&V BP is transformed into a C 3 equivalent by converting the 
new broken-plane C 3flp to an injection velocity at parking 
orbit altitude, adding and converting the sum to a new 

and slightly larger C 3 value at each point on the strip. 

C. Launch/Injection Geometry 

The primary problem in departure trajectory design is to 
match the mission-required outgoing F-infinity vector, V^, to 
the specified launch site location on the rotating Earth. The 
site is defined by its geocentric latitude, <t > 1 . and geographic 
east longitude, X L (Figs. 10 and 11). 

Range safety considerations prohibit overflight of popu¬ 
lated or coastal areas by the ascending launch vehicle. For each 
launch site (e.g.. Kennedy Space Center, Western Test Range, 
or Guiana Space Center), a sector of allowed azimuth firing 
directions Z. is defined (measured in the site's local horizontal 
plane, clockwise from north). For each launch vehicle, the 
allowed sector may be further constrained by other safety 
considerations, such as spent stage impact locations down the 
range and/or down-range significant event tracking capabilities. 

The outgoing K-intinity vector is a slowly varying function 
of departure and arrival date and may be considered constant 
for a given day of launch. It is usually specified by its energy 


magnitude C 3 = |V ( J 2 , called out as C Z L in the plots and 
representing twice the kinetic energy (per kilogram of injected 
mass) which must be matched by launch vehicle capabilities, 
and the K-infinity direction with respect to the inertial Earth 
Mean Equator and equinox of 1950.0 (EME50) coordinate 
system: the declination (i.e., latitude) of the outgoing asymp¬ 
tote 6^ (called DLA), and its right ascension (i.e., equatorial 
east longitude from vernal equinox, T) (or RLA). These 
three quantities are contour-plotted in the handbook data pre¬ 
sented in this volume. 

1. Launch azimuth problem. The first requirement to be 
met by the trajectory analyst is to establish the orientation of 
the ascent trajectory plane (Ref. 7). In its simplest form this 
plane must contain the outgoing V-infinity (DLA, RLA) 
vector, the center of Earth , and the launch site at lift-off 
(Fig. 10). As the launch site partakes in the sidereal rotation of 
the Earth, the continuously changing ascent plane manifests 
itself in a monotonic increase of the launch azimuth. with 
lift-off time, , (or its angular counterpart,- a L , measured 
in the equator plane): 


cotan Z ; 


cos <t>j X tan 5^ - sin 0 L X cos(ft w - o^) 

sin (cr^ - a L ) 
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Fig. 10. Launch/injection trajectory plane geometry 
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Flic dailv tunc lustoiv of a/unuth can be obtained fioin Flq. 
I 1 *) for a given o^, ^ departure asymptote direction by cate- 
full) following quadiant intes explained below and using the 
follow mg auxihaiy expiessions (see Fig 11) 
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A iel.itive launch time, <RI , ■ measuied with lespect to an 
medial leleience (the depailute as\ nipt ole meiidtan's light 
ascension a ), can be defined as 
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Fig. 11. Earth equator plane definition of angles involved 

in the launch problem 

live) and otiO.O deg in the fouith quadiant (when cotan 1 
is negative). 

A generalized plot ot relative launch time f K/ v vs bunch 
azimuth 2J, can be constructed based on Fqs. ( l >) and (12). 
if a fixed launch site latitude is adopted (e.g., 0 7 2N.5 deg 

lot Kenned) Space Flight i'entei at Cape Canaveial. F'lottda). 
Such a plot is piesented in Fig. 12 with dcpuituic asvmptote 
declination as the contoui paiametei. File plot is applicable 
to any icahstic depat lute condition, independent of a v . dale, 
oi Hue launch tune t t (Ret S). 
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lhe expiession lot tlq °) must he used with computa 
titMi.il tegjid tin v|u.»diants. smgulai points, and sign conven- 
lions It the launch is known to be dttect tie. east waul), 
then 1,, when cot.m is negative. must be collected to 
1 1 ,* ( ♦ IS00 1 oi uMiogiade twestwaid) launches. 

IS0 0 deg must be added m the thud (when cotan 2 is post 


2. Daily launch windows Inspect urn of Fig. 12 indicates 
that gcncialt) two contouis exist lot each declination value 
(e.g., JL -10 deg), one occuumg at i h . f r dining the a m. 
hours, the othei m the p.m. bouts of the asvmptote lelatue 

Max" 

♦ 

Since lilt ott tunes aic bounded hv pieseiected launch site 
dependent limiting values ol launch a/limith 2.' ; (eg, "it deg 
and 115 deg), each of the two declination contouis thus con¬ 
tains a segment dining which launch is peinussihlc "a launch 
window.* 4 Fhc two segments on the plot do define the two 
available dailv launch windows. 

\s can be seen liom I ig 12. loi v 0. the two dailv 
launch oppottumties ate sepaiated bv exactlv 12 horns, with 
an meteasmg i thev close m on each othei. until at I 
,0/ - thev hum go into a single dailv oppoituml) Foi |»S J • 
i0/ . a “split" ot that single launch window o.vuis. disallow 
mg an evei mcieasing sectoi ot azimuth values Fins sectoi is 
sviumetuc about east and its hunts can be deteimmed Horn 
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Fig. 12. Generalized relative launch time y vs launch a/imuth and departure asymptote declination b 00m 
Pair of typical example launch windows toi iS^ •* 10 deg are shown by bold cuive segment 

(leproduced from Ref. 8). 



As gets longer. the sector of unavailable launch a/ituuiliN 
reaches the safety houikianes ot permissible launches, and 
planar launch ceases to exist (Fig. U). This subiect \v»ll be 
addressed again in the discussion ot‘ “dogleg'* ascents. 


ft 11*1110 14 in a sketch o! a typical daily launch geometiv 
situation, shown upon a Metcatoi map ot the celestial sphere. 
I he two launch windows exhibit a similar geometry since the 
inclinations o! the ascent tiatectoiv planes are functions ot 
launch site latitude 0 t and a/tmuth only 


cos / - cosV sin (15) 

The two daily opportunities do differ greatly, however, in 
the right ascension of the ast ending mode 11 ot the orbit and 
in the length of the travcised in-plane arc. the range angle f). 


The angular equatorial distance between the ascending 110*’ 1 
and tlu launch site meridian is given by 

sin 6 \ sin I 

sin (a ill “- — (l<>) 

Ouadtjui mles toi tins equation tin oho the observation that a 
negative cos places (a^ - 12> into the second 01 thud quad¬ 
rant, while the sign ot sin (o 7 • 12) determines the choice 
between them. 

I lie range angle f) is measured 111 the menial ascent traiec- 
toiv plane lioni the hft-otf point at launch all the way to the 
departure asvmptote duection. and can be computed lor a 
given launch time i ( 01 - o ! \t i I and an a/imuth L', 

aheady known from Fq ( l M as follows 

cost* - sin\ X sin y* * cos \ \ cos \ costa o ) 
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ig. 13. Permissible regions of azimuth vs asymptote 
declination launch space for Cape Canaveral 


extent of range angle 6 can be anywhere between 
0 deg and 360 deg, so both cos 0 and sin 0 may be desired in 
its determination. The range angle 0 is related to the equa¬ 
torial plane angle, Aa = a m - discussed before. Even 
though the two angles are measured in different planes, they 
both represent the angular distance between launch and depar¬ 
ture asymptote, and hence they traverse the same number of 
quadrants. 

Figure 15 represents a generally applicable plot of central 
range angle 0 vs the departure asymptote declination and 
launch azimuth, computed using Eqs. (17) and (18) and a 
launch site latitude <p L = 28.3 (Cape Canaveral). The twin 
daily launch opportunities are again evident, showing the 
significant difference in available range angle when following a 
vertical, constant 6 M line. 

It is sometimes convenient to reverse the computational 
procedure and determine launch azimuth from known range 
angle 0 and t RlT , i.e., Aa = - a L , as follows: 
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cos 6 m X sin (a. ~ o^) 
sin 6 


cos £, = 


sin - cos 6 X sin 
sin 0 X cos <t> L 


with the lighting conditions at lift-off and consequently allows 

(19) a lighting profile analysis along the entire ascent arc. 

The angle ZALS, displayed in Fig. 16, is defined as the 
angle between the departure V M vector and the Sun-to-Earth 

(20) direction vector It allows some judgment on available ascent 
lighting. 


Figure 16 displays a 3-dimensional spatial view of the same 
typical launch geometry example shown previously in map 
format in Fig. 14. The difference in available range angles as 
well as orientation of the trajectory planes for the two daily 
launch opportunities clearly stands out. In addition, the figure 
illustrates the relationship between the “first” and “second 
daily” launch windows, defined in asymptote-relative time, 
r RiT' as ^ntrasted with “morning” or “night” launches, 
defined in launch-site-local solar time. The latter is associated 


The length of the range angle required exhibits a complex 
behavior- the first launch window of the example in Figs. 14 
and 16 offers a longer range angle than the second, but the 
second launch window opens up with a range angle so short 
that direct ascent into orbit is barely possible. Further launch 
delay shortens the range even further, forcing the acceptance 
of a very long coast (one full additional revolution in parking 
orbit) before transplanetary departure injection. A detailed 
analysis of required arc lengths for the various sub-arcs of the 
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Fig. 18. Angle from perigee to departure asymptote 


v { = is the true anomaly of the injection point, usually 
is near 0 deg. and can be computed by iteration 
using: 


tan 7) 


t’ H sin 

1 + e H cos v f 



y g = injection flight path angle above local horizontal, 
deg 

v = eccentricity of the departure hyperbola: 



1 + 





(24) 


To make Eq. (22) balance, the parking orbit coast arc, 
must pick up any slack remaining, as shown in 
Fig. 17. A negative 0 coaST implies that the -solution was 
too short-ranged. A direct ascent with positive injection true 
anomaly v g tie., upward climbing flight path angle, 7 /f at 
injection) with attendant sizable gravity losses, may be accept¬ 
able. or even desirable (within limits) for such missions. 
Alternately, the other solution for . exhibiting the longer 
range angle 0 , and thus a longer parking orbit coast. ^coast* 
should he implemented An extra ievolution in parking orbit 
may be a viable alternative. Other considerations, such as 
desire for a lightside launch and/oi injection, tracking ship 
location and booster impact constraints, may all play a signif¬ 


icant role in the ascent orbit selection. A limit on maximum 
coast duration allowed (fuel noil-off, battery life, guidance 
gym drift, etc.) may also influence the long/short parking orbit 
decision. In principle, any number of additional parking orbit 
revolutions is permissible. Shuttle launches of interplanetary 
missions (e.g., Galileo) are in fact required to use such addi¬ 
tional orbits for cargo bay door opening and payload deploy¬ 
ment sequences. In such cases, however, the precessional 
effects of Earth’s oblateness upon the parking orbit, primarily 
the regression of the orbital plane, must be considered. 

4. Parking orbit regression. The average regression of the 
nodes (i.e.. the points of spacecraft passage through the 
equator plane) of a typical direct (prograde) circular parking 
orbit of 28.3-deg inclination with the Earth's equator, due to 
Earth's oblateness, amounts to about 0.46 deg of westward 
nodal motion per revolution and can be approximately com¬ 
puted from 



540° X r 2 X / X cos i 

f 2 



() 


deg/revolution (25) 


where 

r ~ Earth equatorial surface radius. 6378 km 
r - circular orbit radius, typically 6748 km for an 

O ■ 

orbital altitude of 370 km (200 nmi) 
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J 2 - 0.00108263 for Ear lit 

i = parking orbit inclination. dc i» * computed for 
a given launch geometry (v y 



(26) 


This. cuirectiAn. 
revolutions, must In 
time and, hence, t 


multiple* 1 ■/ .he orbital stay time of A 
; ronsid. »» 2 it; determining a biased launch 
he ngh* ;jtension of the launch site at 


lift-off: 


o. - u, + il X A f deg (27) 

eW 1 ‘ 

5. Dogleg ascent. Planar ascent has been considered exclu¬ 
sively, thus far. Reasons for performing a gradual powered 
plane change maneuver during ascent may be many. Inability 
to launch in a required azimuth direction because of launch 
site constraints is the prime reason for desiring a dogleg ascent 
profile. Other reasons may have to do with bum strategies 
or intercept of an existing orbiter by the ascending spacecraft, 
especially if its inclination is less than the latitude of the 
launch site. Doglegs are usually accomplished by a sequence of 
out-of-plane yaw turns during first- and second-stage bum, 
optimized to minimize performance loss and comment* ig as 
soon as possible after the early. low-altitude, high aerodynamic 
pressure phase of flight is completed, or after the necessary 
lateral range angle offset has been achieved. 

By contrast, powered plane change maneuveis out of 
parking orbit or during transplanetary injection are much less 
efficient, as a much higher velocity vector must now be rotated 
through the same angle, but they may on occasion be opera¬ 
tionally preferable. 

As already discussed, a special geometric situation develops 
whenever the departure asymptote decimation magnitude 
exceeds the latitude ot the launch site, causing a "split azi¬ 
muth" daily launch window Figure 13 shows the effects of 
asymptote decimation and range safety constraints upon the 
launch problem As the absolute value ot declination increases, 
it eventually reaches the safety constraint on azimuth, prevent¬ 
ing any further planar launches The situation occurs mostly 
early and or late in the mission's departure launch period, and 
is frequently associated with dual launches, when month-long 
departuie periods are desired 

The Shuttle-era Space Transportation System (STS), includ¬ 
ing contemplated upper stages, is capable of executing dogleg 
operations as well. These would, however, effectively reduce 
the launch vehicle's payload (or C\l capability, as they did on 
expendable launch vehicles ot the past. 


6. Tracking and orientation. As the spacecraft moves away 
from the Earth along the asymptote, it is seen at a nearly 
constant declination that of the departure asymptote, 6 M 
(DLA in the plotted data). The valrc ot D! A greatly affects 
tracking coverage by stations located at various latitudes; 
highest daily spacecraft elevations, and thus best reception, are 
enjoyed by stations whose latitude is closest to DLA. Orbit 
determination, using radio doppler data, is adversely affected 
by DLA's near zero degrees. 

The spacecraft orientation in the first few weeks is often 
determined by a compromise between communication (antenna 
pointing) and solar heating constraints. The Sun-spacecraft- 
Earth (SPE) angle, defined as the angle between the outgoing 
K-infinity vector and the Sun-to-Earth direction, is very 
useful and is presented in the plots under the acronym ZALS. 
It was defined in the discussion of Fig. 16. This quantity has 
many uses: 

(I) If the spacecraft is Sun-oriented, ZALS equals the 
Earth cone angle (C'A), depicted in Fig. 19 (the cone 



Fig. 19. Definition of cone and dock angle 
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an jngle between ilic vehicle s longitudinal --/a\is..nd 
the object direction). 

(2) The Sun-plase angle. <b v (phase angle is the Sun-object - 
spacecraft angle) describes the state of the object's 
disk lighting: a fully lit disk is at zero phase. For the 
Farth (and Moon), several days after launch, the su.i- 
phase angle is: 

< 1 * - IMI-/AIS [IS) 

4 ■ 

(3) The contour labeled ZALS ” 4 H)'* separates two cate¬ 
gories of tianster trajectoiies those eail> depaituies 
that first cut inside Faith's mbit, thus stalling out at 

negative heliocentuc true anomahe. toi ZALS >'>i) 

% 

and those latei one* that stait at positive tiue anoma¬ 
lies. heading out toward Jupitei and never exper tcik mg 
the tncieased soldi heating at distaiues of less than 

1 At tor ZALS v. **0' 


Annliritinfu In.* I^Atl cfii.k ( R»»t U| r.iititw fn nl’ 
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the techniques available, such as passive wait tor natural 
alignment of the continuously regressing .space station orbit 
plane (driven by Farth’s oblateness) witn the required 
’’•infinity vector, or the utilization of 2- and 3- impulse man¬ 
euvers. seeking to perform spacecraft plane changes near the 
apogee of a phasing orbit where velocity is lowest and thus 
turning the orbit is easiest. These two approaches can be 
combined with each other, as well as wi»h other suitable 
maneuvers, such as. 

(I) IVep space ptopulsive burns tor orbit shaping and 
phasing. 

il) (iruvitv assist tlvbv. includim* Farth return AYFXiA, 

131 Aoiody iij*iic turns at grazing perigees oi at inter¬ 
mediate planetary s»» ingbys.and 

(4) Multiple revolution intention burns, requiring several 
low. gia/ing passes, combined with apogee plane change 
maneuvers, etc. 


7. Post-launch spacecraft state. Alter the spaceman has 
departed fiom the immediate vicinity ot Faith li.c.. lett the 
laith's qiheie ot influence ■*! about 1-2 million km), it moves 
on a heliocentric conic whose mtii.il conditions may be 
approximated as 



V 


/ is in 


f V 


l>M 
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f V 


.V C 


A/ 


(*U) 


wlieie R and V ot i jith aie evaluated iron 1 an ephcmeiN at 
tune ot miectum and A/ lepiescuts time elapsed Miue then (in 
seconds) 1 lie V vovlor m I Ml 50 vattcsian vooidmale van 
be vMivtiucteJ usiui* \ (\i \ Dl A - S and Kl A o 

ill tlnee k.o|oi*«»iietits as follows 
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H. Orbital launch problem Orbital launch tiom the Shuttle, 
ttom other elements ot the SIS -u trom ativ tempoiarv or 
permanent oihital space station complexes, mtioduces eiittielv 
new concepts uito die I arth-departure ptohlem Some ot 
the new k''listiamts. .ilteadv mentioned, limit ou: abihtv 
to launch a given uiteip. metaiv mission llvsiowlv tegiessing 
space station othn (see Fq 2>) genet ally docs not contain 
tlie I intimtv ve^ioi requited at depaitme Oibit litetime 
oi othei soibideiatlons ma> dktate a .pave station's orbital 
altitude that mav be too high tor an elticienl intention bun 
Innovative depaiture strategies jie bee m nine to erne tee. 
attempting to alleviate these problems a ievent Svieine 


All o| these devices can be optimized to permit satisfactory 
orbital launches, as well as to acliieve the most desirable condi¬ 
tions at the final anival body. In general, space launch advan¬ 
tages. such as ou-oibit assembly and chevkout of payloads and 
chisVrcd multiple propulsion stages, or oihital const ruction of 
bulky and t-agile subsystems (solar panels, sails, aiiiennas. 
radiatois. booms, i.c.l will, it is hoped, greatly outweigh the 
significant deep-space mission penalties incurred because o! 
the space slat mil's inherent oihital orientation mconipalihilitv 
with depaiture requirement 


D. Planetary Arrival Synthesis 

I lie planctaiv arrival traject »i\ design problem involves 
satistving the project s cngmccung and science objectives at 
lire target body by shaping the amval tiajectoiv m a citable 
mannei. As these objectives mav be quite diverse. only lo*:t 
illuslialive svcnanos shall be discussed i.i fills section flv 1 v . 
oibitei atumspheiiv piohe. and to a small extent .ai.dei 
missions 


I Flvbv trajectory design. In this mission mode, be amval 
tiajectoiv is not modified m am deterministic wav at the 
planet the oiigm.il aim point and amval time are vhosen to 
satistv trie largest ntmibei ot potential objectives. long before* 
band 

Ibis pto„ess involves the choke i l arrival dak to ensuie 
desirable characteristics. siu'i j\ i l »e values ot the variables 
Mil*. \*\\* /APS. ck . presented in plotted toim m (he data 
(ion ot tins \olume 
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DAP, (lie planet I'tju.it'Mi.il iUvlm.ition. <S_ . ot tin* incoming 
asunptole. uv. ol the 1 intiniu vectoi. pi ovules ihe measine 
ot ilu* minimum possible inclination ot ll\b\ Its negative is 
also known as the latitude ot \ettieal uup.u I (1 VI). 

The magnitude ol I intiiiit\, VI IP jY J. enables one to 
couliol tlu* tl>b> turn angle AC between the incoming atul 
outgoing V m veetois b\ a suitable elioiee ol elosest appioaeh 
(<' *U i.ulms,/ (see lit* iU 

AC' - ISO -V . deg < i.') 

wheie the asNinptote halt angle, is Iouiul horn 


eos/i % - (Ml 

' I * t 

. . "* /• 


VIII* also enables the destgnet to evaluate planetoeentue 
velocity, i . at anv distance, /. on the Ilyin hypeibola 

//i t 

I I J * * 1 . km s 

In the e.hcvc equations. jj ;i (oi (».W ). in the giavitational paiam- 
etei o* the attival bod\ , 
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Another pair of significant variables on which to base 
arrival date selection are ZAPS and ZAPE-the angles between 
V infinity and the planet-to-Sun and -Earth vector* respec¬ 
tively. These two angles represent the cone angle (CA) of the 
planet during the far-encounter phase for a Sun- or an Earth- 
oriented spacecraft, in that order. ZAPS also determines the 
phase angle. 4> s , of the planet's solar illumination, as seen by 
the spacecraft on its far-encounter approach leg to the planet: 

* s - 180-ZAPS (35) 

Both the cone angle and the phase angle have already been 
defined and discussed in the Earth departure section above. 


The flyby itself is specified by the aim point chosen upon 
the arrival planet target plane. This plane, often referred to 
as the 5-plane, is a highly useful aim point design tool. It is a 
plane passed through the center of a celestial body normal to 
Y^, the relative spacecraft incoming velocity vector at infinity. 
The incoming asymptote, i.e.. the straight line, zero-gravity 
extension of the V . -vector, penetrates the 5-planc at the aim 

point. This poinl.de.d by the target vector B in the 5-plane, 

is often described by its two components B * T and B * R, 

A A 

where the axes T and R form an orthogonal set with V^. The 
T axis is chosen to be parallel to a fundamental plane, usually 
the ecliptic (Fig. 2i) or alternatively, the planet's equator. The 
magnitude of B equals the semi-minor axis of the flyb\ hyper¬ 
bola. b. and can be related to the closest approach distance, 
also referred to as the periapse radius, r , by 


or 



oo 



^ ■ r ■ 


km (36) 



1/2 


+ IBI 


. km (37) 



B MISS PARAMETER. BIS 
(TARGET VECTOR) 

6 AIM POINT ORIENTATION 
S PARALLEL TO INCOMING ASYMPTOTE. V M 

T PARALLEL TO ECLIPTIC PLANE 
AND 1 TO S' 

R =Sx T 


Fig. 21. Definition of target or arrival 5-plan* coordinates 


The two systems of 5-plane T-axis definition can be recon¬ 
ciled by a planar rotation. -A0. between the ecliptic T- and 

A A, A, 

R-axes and the T- and R-axis orientations of the equator 
hased system 


The direction angle 0 of the 5-vector, B. measured in the 
target plane clockwise from the T axis to the 5-vector position 
can easily be related to the inclination, i, of the flyby trajec¬ 
tory. provided that both (DAP) and the T-axis, from which 
0 is measured clockwise, are defined with respect to the same 
fundamental plane to which ‘he inclination is desired. For a 
system based on the planet equator (Fig. 22): 

cos i.,... = cos 0 V , X cos (38) 

try PK q 

which assumes that Q PF q is computed with the T axis parallel 
to the planet equator (i.e.. Tp F Q = X POLEp^k not the 
ecliptic, as is trequentK assumed (T^ t7 = X POLE^j). 
Care must be taken to use the 0-augle as defined and intended. 


sin (a. - a h:p ) 

tan A0 -- 1 ———;- . - -:-: 

cos 6 X tan o,.„ - sin o X cos (a - ot rD ) 

°" r.r °® t.r 

(39) 

where 

a^p and h Ep are right ascension and declination of the 
ecliptic pole in planet equatorial coordinates. For Jupiter, 
using constants in Section V: 

a Fp * 290.598 b Ep = 87.789, deg 

<K. and h m are RAP and DAP, the directions of incoming 
V „. also in planet equatorial coordinates. 
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ECLIPTIC POLE IEPI 



« PLANE (1 TO V.J 


T / ( 7 " V « X P H I. 

t /‘AU * v<x ' x p 7 ’/V 

Aa--90- I {a tp - 90) - (a*, - 180)1 - a f7> 

AO ARCTAN (SIN Aiv (COS ^ x TAN f> Fp SIN fi x , x COb Aft)l 

Fig. 22. Two T axis definitions in the arrival #-plane 


b 

k 


The correction AO is applied to a 0 angle computed in the 
ecliptic system as follows (l\e ”1: 

( V/ y» ' <',ct * At ' (401 

The ecliptic T. K rcl axes, how cum. have to be totaled b> 
(clockwise diiechon is positive in the tfplane) to obtain 
planet equatonal T. cooidmate axes The ^-magnitude 

of an aim point in eithei s\ stem is the same 


pei iapse, at the entiv point ot a pmbe, 01 genet ally at an> 
position r (subsetipt .V --- Sun. could be teplaced by l ~ Farth. 
if desired), (see Fig. .'4): 


cos *1* “ 

4 * 


-cos fL \ cos / Al’ v - sin ^ \ sin / Al’ v 


\ cos (/•/'/’-0 ) 


(41) 


vvheie 


The pio|cction$ of the SuiMo-planct and Faith-to-planet 
vectors into the H plane represent aim point loci of diaiuetiic 
Sun and Faith ivcul tat ions, lespectivelv, as defined in Fig. 25. 
The H plane angles (w ith tespect to T fVi axis) of these van* 
a hies ate presented and labeled FTSP and FTFP, respectivelv, 
in the plotted mission data. In addition to helping design ot 
else avoid dianietnc occupations. these i|iiantittes allow 
computation *4 phase angles. «h fl of the planet at the spaceciaft 


c• r 1,10 Jl,n P omt anglej-.i the H plane, must be w ith 
icspect to the same T. as FTSP 

^ - the amval unge ingle fioin tntinitx (a position tai 
out on the incoming asvmptotc) to the point ot 
mteiest r (Fig 25) 

I he computation of the auival »ange angle, to the 
position ot the desited event depends on its t\pe. as follows 
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ARRIVAL 
/? PLANE 


SUN (OR EARTH) TO 
PLANET VECTOR, 


SHADOW 
CYLINDER 
(SUN OR 
EARTH) 



TO SUN 
(OR EARTH) 


I. At flyby periapse. (/?*, > VO deg): 


= cos(-i^) 




(r p is periapst radius). 

2. At given radius r . anywhere on the flyby trajectory 
(see Fig. 25): 


ft. = + ", H3) 

^ should be computed from periapse equation. 
Eq. (42) v f at r can be obtained from (- f b < v f < + 
v f has negative values on the incoming branch): 




At the entry point hiving a specified flight path angle 
(Fig. 251. 


Fig. 23. Definition of approach orientational coordinates 
ZAPS and ETSP. ZAPF. and ETEP 
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Fig. 24. Phase angle geometry at arrival planet 
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Fig. 25. Typical entry and flyby trajectory geometry 
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where v F is the true anomaly at entry, should always 
he negative, and can be computed if entry radius and 
altitude, r E = r SVRF + h F and y t:% the entry angle, are 
known: 



whereas the fictitious periapse radius r p for the entry, 
to satisfy y E at r t> is equal to 



The B value corresponding to this entry point can be com¬ 
puted from Eq. (36), while t l, e 0 .. de in the £-plane would 
depend on the desired ent*., latitude inclination, Eq. (38), or 
phase angle, i.q. (41). 

The general flyby problem poses the least stringent con¬ 
straints on a planetary encounter mission, thus allowing 
optimization choices from a l.nge its: %.f secondary parameters, 
such as satellite viewing and occultation. planetary fieUs and 
particle in situ measurements, special phase-angle effects, etc. 

A review of the plotted handbook variables, required in the 
phase-angle equation (Eq. 41), shows that the greatest magni¬ 
tude variations are experienced by the ZAPS angle, which is 
strongly flight-time dependent: the longer the trip, the smaller 
ZAPS. For low equatorial inclination, direct flyby orbits, this 
implies a steady move of the periapse towards the lit side and. 
eventually, to nearly subsolar periapsesfor long missions. This 
also implies that on such flights the approach legs of the tra¬ 
jectory are facing the morning terminator or even the dark 

side, as trip time becomes longer, exhibiting large phase angles 
(recall that phase is the supplement of the ZAPS angle on 
the approach leg). This important variation is caused by a grad¬ 
ual shift of the incoming approach direction, as flight time in¬ 
creases. from the subsolar part of the Jovian leading hemispheic 
(in the sense of its orbital motion) to its antisolar part. 


The arrival time choice on a very tine scale may greatly 
depend on the desire to observe specific atmospheric/surface 
featuics (e.g., the Great Red Spot) or to achieve close encoun¬ 
ters with specific satellites of the arrival planet. Passage* 
through special satellite event zones, e.g.. flux tubes, wakes, 
ge«ventric and or heliocentric occupations, require close 


control of arrival time. The number of satellites passed at 
various distances also depends on the time of planet CIA, 
These fine adjustments do, however, demand arrival time 
accuracies substantially in excess of those provided by the 
computational algorithm used in this effort, which generated 
the subject data (accuracies of 1-5 min for events or 1-2 h 
for encounters would be reo 1 ired vs uncertainties of up to 
1.5 days actually obtained with the rectilinear impact pseudo¬ 
state theorem). Numerically searched-in integrated trajec¬ 
tories. based on the information presented as a first guess 
input, are mandatory for such precision trajectory work. 

Preliminary design considerations for penetrating, grazing, 
or avoiding a host of planet-centered fields and particle struc¬ 
tures. such as magnetic fields, radiation belts, plasma tori, ling 
and debris structures, occupations by Sun. Earth, stars, or 
satellites, etc., can all be presented on specialized plots, e.g.. 
the £-plane. and do affect the choice of suitable aim point 
and arrival time. All of these studies require the propagation 
of a number of flyby trajectories. Adequate initial conditions 
for such efforts can be found in the handbook as: VHP ( V ao ) 
and DAP (6^) already defined, as well as RAP (a^). the 
planet eqi atoria' right ascension of the incoming asymptote 
(i.e.. its east longitude from the ascending node of the planet's 
mean orbital plane on its mean equator, both of date). The 
designer's choice of the aim point vector, either as B and 8 , or 
as cartesian B * T and B * R. complete the input set. Suitable 
programs generally exist to process this information. 


2. Capture orbit design. The capture problem usually 
involves the task of determining what kind of spacecraft orbit 
is most desired and the interconnected problem of how and at 
what cost su:h an orbit may be achieved. A scale of varying 
complexity may be associated with the effort envisioned an 
elliptical long period orbit with no specific orientation at the 
trivial end of the scale, through orbits of controlled or opti¬ 
mized lines of apsides (i.e.. periapse location), nodes, inclina¬ 
tion. or a safe perturbed orbital altitude. Satellite G/A-aided 
capture, followed by a satellite tour, involving multiple satel¬ 
lite G/A encounters on a numher of revolutions, each designed 
to achieve specific goals, probably rates as the most complex 
capture orbit class Some orbits arc energetically very difficult 
to achic/e, such as close circular orbits, but all require signifi¬ 
cant expenditures of fuel. As maneuvers form the background 
to this subject a number of useful orbit design concepts shall 
be presented to enable even an unprepared user to experiment 
with the data presented. 

The simplest and most efficient mode of orbit inaction is 
a coplanar bum at a common periapse of the arrival hyperbola 

and the resulting capture orbit (Fig. 26). The maneuver AF 
required is: 
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2 H \ r t 
r [r i r ) 

r ' I /' 


, km n (4S) 


A plot ot' orbit inset (ion AI u\| lined as a function of r /t ami 
/* iusing r*|. 50) is piesented in l ; ig. ?\ Hie apoapse ladius of 
such an oibit ot given peiuul would be 


The oibital peiuul foi such an oibit. tcquiiing knowledge 
of penapse and apoapse ladii. r t and . is 



It’ on the othei hand, a known oibit peiuul /* (in seconds) 
isdesued. the exptession tot AI is 



/ , 

/:/i t \ /■• 

r, - x/ : — - - r r .kin Pli 

An evaluation ot* l\|. (50) (and lig. -shows that lowest 
oibit niseition AI is obtained t’oi the lowest value of r the 
longe,t peiuul 1\ and the lowest l \. ot auival 

Ot some inteiest is imcction into citcular captuie mbits, a 
special case of the coapsulal insertion pioblem. It can be 
shown (Ret. S) that an optimal A I* exists foi inset lion into 
captuie oibits of constant eccentiicitv. including »* - 0. u\. 
cuculai oibits. which would icquitc a specific ladius 



Fig. 26. Coapsidal and cotangentnl capture orbit insertion geometries 
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Fig. 27. Coapsidal capture orbit insertion maneuver AV 
requirements for Jupiter (using Eq. 50) 


Aco p , one can solve for the hyperbolic periapse r CA and the 
burn radius r using selected values of true anomaly at hyper¬ 
bolic burn point v t1 and !♦« capture orbit equivalent 


p.. = v • + Ac o 

f- it p 


(55) 


utilizing the following three equations (where E and H stand 
for elliptic and hyperbolic, respectively): 



and 





while the corresponding optimal value for A V would be 



! 




Frequently the orbital radius obtained by use of Eq, (52) is 
incompatible with practical injection aspects or with arrival 
planet science and engineering objectives. 


A more general coplanar mode of capture orbit insertion, 
requiring only tangentiality of the two trajectories at an 
arbitiai\ nun-, aver point of radius r common to both oibits. 
Fig 20. requires a propulsive effott of 


Al 




»n {r A 


r(r 


A 


+ r - r) 

p 


(54) 


It can be clearlv seen that h> performing the burn at periapse 
the substitution r = brings us back to Eq. (48). 

The cotangential maneuver mode provides nonoptimal con¬ 
trol mot the orientation of the major axis of the capture orbit. 
If it i\ desired to rotate this line of apsides clockwise by 


The procedure of obtaining a solution to these equations is 
iterative. For a set of given values for r ,, r . and an assumed 

* p 

Aco. a set of \' H and i^. the hyperbolic and elliptical bum 
point true anomalies which would satisfy Eqs. (55-58) can 
be found. This in turn leads to r, the maneuver point radial 
distance, and hence. AT (Eq. 54). A plot of A I'cost for a set 
of consecutive Aco^ choices will provide the lowest AC value 
for this maneuver mode. 

For an optimal insertion into an mbit of an uihiii ji\ niajoi 
axis orientation one must turn to the moie general, still co¬ 
planar. but intersecting (i.e., noutangeniial bum point) man- 
euvei (see Fig. 28). It provides sufficient flexibilit) to allow 
numerical optimization of AT with respect to apsidal rotation. 
A co - 

p 

A more appropriate wa> to define apsidal orientation is to 
measure the post-maneuver capture orbit periapse position 
angle with respect to a fixed direction, e.g.. a far encounter 
point on the incoming asymptote. -V M . thus defining a cap¬ 
ture oihit periapse range angle 
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where 


HYPERBOLIC FLYBY 
PERIAPSE,r C/1 .CAN 

VAR Y: r A >r CA >0, 

PROVIDED r $ IS AVOIDED. 

OPTIMAL BURN AV 

REQUIRES t nA >r 
C A p 


geometries 


- Ac J 


+ arc cos /--\ 


(59) 


Taken trom Ref. 8. the expression for the intersecting burn 
•il'is 

ai j . I* |i- —L_]. J. r*T V/> 

71 L r %J L 



) 

( 60 ) 


r ~ planet-centered radius at hum. r ->r • r 

*■ (.4 

r x and r f . = anoapse and penapse radii ol capture ellipse 


n 


* 
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r cA = closest approach radius of flyby hyperbola 

$ * flag: 6 = +1 if injection occurs on same left 

(inbound or outbound) of both hyperbola 
and capture ellipse, 6 = -1 if not. 


It should be pointed out that Eq. (57) and (58) still apply 
in the intersecting insertion case, while Eq. (56) does not 
(as it assumes orbit tangency at burn point). 

The evaluation of intersecting orbit insertion is more 
straightforward than it was for the cotangential case. By assum¬ 
ing K. an d orbit size (e.g., F = 6 km/s, r = 4 RJ, P - 200 

P 

days, similar to a Galileo capture orbit) and stepping through a 
set of values for v £ . the capture orbit burn point true anomaly, 
one obtains, using Eqs. (57-60), a family of curves, one for 
each value of R CA , the hyperbolic closest approach distance. 
As shown ip. Fig. 20, the envelope of these curves provides the 
optimal insertion burn AT for any value of apsidal rotation 
Ato,, desired. The plot also shows clearly that cotangential and 
apsidal insertion burns arc energetically inferior to burns on 
the envelope locus. For the same assumed capture orbit, a 
family of optimal insertion envelopes, for a range of values 
of arrival F^. is presented in Fig. 30. The smallest value. 
V m = 4.458 km/s, represents the remaining pre-insertion 
energy after an lo encounter, typical of Galileo mission 
strategy. 


The location of periapse with respect to the subsolar point 
is of extreme importance to many mission objectives. It can be 
controlled by choice of departure and arrival dates, by AF 
expenditure at capture orbit insertion, by an aerodynamic 
maneuver during aerobraking, by depending on the planet's 
motion around the Sun to move the suhsolar point in a manner 
optimizing orbital science, or by using natural perturbations 
and making a judicious choice of orbit size, equatorial inclina¬ 
tion, r, and initial argument of periapsis. w n , such as to cause 
regression of the node, 12, an J the advance of periapsis, cb. 
both due to oblateness to mo/e the orbit in a desired manner 
or at a specific rate. Maintenance of Sun-synchronism could 
provide constant lighting phase angle at periapse, etc., by some 
or all of these techniques. For an elliptical capture orbit 
(Fig. 31). 


il 



cos I X 



( 61 ) 


U? 




- (5/2) sin 1 


fix 


180 

rr 


. deg/s 



n 



, mean orbital motion, rad/s 


(63) 


A P 

a = —- — , semi-major axis of ellipitical orbit, km (64) 

T P 

p = - - — — , semi-latus rectum of elliptical orbit, km (65) 

A P 

R s = Jupiter's equatorial surface radius, km 

J 2 = Jupiter's oblateness coefficient (for values see 
Section V on constants). 

For the example orbit (4 a 272.4 RJ) used in Figs. 28 
and 29, if near equatorial, the regression of the node and 
advance of periapse, computed using Eqs. (61-65), would 
amount to a negligible -0.23 and +0.46 deg/year, respectively. 
For contrast, for a grazing (a = 71,398 km) low-inclination, 
near-circular orbit, the regression of the node would race along 
at -64.5 deg/day, while periapse would advance at 129 d?g/ 
day, i.e., it would take less than 3 days for the line of apsides 
to do a complete turn. 

It should be noted that il * 0 occurs for i = 90 deg, while 
(1) = 0 is found for i = 63.435 deg. Sun-synchronism of the 
node is achieved by retrograde polar orbits, e.g., a 5 RJ cir¬ 
cular orbit should be inclined 111.08 deg. 

3. Entry probe and lander trajectory design. Entry tra¬ 
jectory design is on one hand concerned with maintenance of 
acceptable probe entry angles and low relative velocity with 
respect to the rotating atmosphere. On the other hand, the 
geometric relationship of entry point, subsolar point and 
Earth (or relay spacecraft) is of paramount importance. 

Lighting during entry and descent is often considered the 
primary problem to he resolved. As detailed in the flyby and 
orbital sections above, the choice of trip time affects the value 
of the ZAPS angle which in turn moves the entry point tor 
longer missions closer to the subsolar point and even beyond, 
towards the morning terminator. 

Landers or balloons, regardless of deceleration mode, 
prefer the morning terminator entry point which provides a 
better chance for vapor-humidity experiments, and allows a 
longer daylight interval for operations following arrival. 

The radio-link problem, allowing data flow directly to 
Earth, or via another spacecraft in a relay role, is very complex. 
It could require studies of the Earth phase angle at the entry 
locations or alternately, it could require detailed parametric 


where 
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Fig. 29. Characteristics of intersecting capture orbit insertion and construction of optimal bum envelope at Jupiter 
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Fig. 30. Minimum SV required for insertion into Jupiier capture orbit of given apsidal orientation 
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Fig. 31. General satellite orbit parameters and precessional motion 
due to oblateness coefficient J 2 (from Ref. 8) 


studies involving relative motions of probe and relax spaceciatt 
throughout probe entry and its following slow descent. 

Balloon missions could also involve consideration of a 
variety of wind dntt models, and thus, ate even more complex 
a., tar as the communications problem with the Farth oi the 
spacecraft is concerned 

E. Launch Strategy Construction 

1 he oMisuaints and desires, briefly discussed above, max be 
displayed on the mission space launch arnxal day plot as 
being limited by the voiitom boundaries «»l the dates 

1)1 A. VHP. ZAP eu thus displaying the allowable lauiuli 
spa^e 

Within this laurivh space a preferred dax by day launJi 
strategx must be specified in accordance with prevailing 
objectives I he simplest lauruh strategx .otten used to maintain 
a constant anixal date at the '.aieet planet results in dailv 
launch points on a fion/onta! line li**m leftmost to rightmost 
maximum allowable ( J bouriuarx tor that aiiivjl date 
Such a strateex makes use ot the tact that most anixal chaiac 
(eristics max >rax nearlx constant across the launch space 
Lighting and satellite positions in this case are fixed, thus 
allow me a similar encounter satellite t» V or satellite tour 


A different choice ot strategy could be to follow a contour 

lin«* lit ftU . • li-i r• tii-'li -it 1)1 \ * it / *«%••! I 

MMC % « It MIC ^ Ot»« t«s « w 4 • n « l«« V t « do I * | .1 I 4 0 « M kl , ^ iflX Sft 


also follow the minimum value locus ol a parameter, e.g.. 
(y. (t.e . the boundary between C lass I and within Type I 
or II) tor each launch dale, throughout the launch space. 


I undamciiUlly different is a launch strategy tor a dual oi 
multiple spacecratt mission, involving more than one launch, 
either ot which may possibly pursue divergent objectives As 
an example Fig M shows the Voyagei. 1 and 2 launch strut 
egy, plotted on an Farth departure vs Saturn arrival date plot. 
A 14-day pad turnaround so pare' n between Punches was 
to be maintaiiied. a 10-day opportunity was to be available 
for each launch, and the two spacecraft had substantially 
different objectives at Jupiter and Saturn one was to be 
lo-mtensive and a close J.ipitci flyby, to be followed by a 
close Titan encounter at Saturn, and the other was Ganymcdc- 
and or Cal listo-in tensive. a dist •::» Jupiter flyby, as a safety 
precaution agams: Jovian radiation damage, aimed to continue 
past Saturn 10 Uranus and Neptune. Here, even the spacecraft 
departure order was reversed bx the strategy within the 
launch space. 


Launch strategies for oibita! departures from a space station 
in a specific orbit promise to introuvice new dimensions into 
mission planning and design. New concepts are beginning to 
emerge on this subject e.g.. Refs. 9 and 10. 


o • ■. • »«*'• * •*. * 



Fig. 32. Voyager (MJS77) trajectory space and launch strategy 
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IV. Description of Trajectory Characteristics 
Data 

A* General 

The data represent trajectory peiformance information 
plotted in the departure date vs arrival date space, thus defin¬ 
ing all possible direct ballistic transfer trajectories between the 
two bodies within the time span considered for each oppor¬ 
tunity. Twelve individual parameters are contour-plotted. The 
first, C 3 /,, is plotted bold on a Time of Flight (TFL) back¬ 
ground; the remaining ten variables are plotted with hold con¬ 
touring on a faint C 3 /. background. F.leven plots arc presented 
for each of twenty mission opportunities between I9«5 and 
2005. 

The individual plots are labeled in the upper outei corner 
by bold logos displaying an acronym of the variable plotted, 
the mission's departure year, and a symbol of the target planet. 
These permit a quick and fail-safe location of desired informa¬ 
tion. 

B. Definition of Departure Varabtes 

C,/.: Farth departure energy (km 2 /s 2 ); same as the 
square of dt^ allure hyperbolic excess velocity 
l f ’ 3 L = V) - 2 %//?,, where 

V g = conic injection velocity (km,'s). 

R t = R s + injection radius (km), sum of 
surface radius R Spj Y/ . f . and injection 
altitude hg , where Rs hA R rf/ refers to 
Earth’s surface radius. (For valu*, see 
Section V on constants.) 

^ 38 gravitational constant times mass of # he 
launch body (for values, refer to Section V 
on constants) 

CjL must be equal to or exceeded by the launch 
vehicle capabilities. 

DU: h m . l . geocentric declination (vs mean Farth equator 
of 1950.0) of the departure \\ vector. May im¬ 
pose launch constraints (deg). 

RLA: , geocentric right ascension (v\ mean Farth 

equator and equinox of 1 1 >S().0) °f the departure 
I' vector. Can be used with C\i. and DLA to 
compute a heliocentric initial state for traiectory 
analysis (deg) 

ZALS: Angle between depature vector and Sun-Faith 
vector. Equivalent to Farth-prohe-Sun angle several 
days out (deg). 

C. Definition of Arrival Variables 

V iIP l * . planetocentnc arrival hvnerbolu excess 
veloc:t> or l-infinity (km s). the magnitude of the 


vector ohtained hy vectoiial subtraction of the 
heliocentric planetary orbital velocity from the 
spacecraft arrival heliocentric velocity. It repre¬ 
sents planet-relative velocity a* -reat distance from 
target planet, at heginning o, fur encounter. Can 
he used to compute spacecraft velocity at any 
point r of flyby, including CIA (periapse) dis¬ 
tance r p : 



where 

Vpfjrprrnk ~ R rav * tationa * parameter GM 
sys ppm \ of the arrival planet system- 
Jupiter plus all satellites. (For 
values, refer to Section V on 
constants.) 

DAP planetocentric declination (vs mean plane* 

equator of date) of arrival vector. Defines 
lowest possible flyby/orbiter equatorial inclination 
(deg). 

RAP a* 4 , planetocentnc right ascension (vs mean 
planet equator and equinox of date, i.e., RAP is 
measured in the planet equator plane from ascend¬ 
ing node of the planet's mean orbit plane on the 
planetary equator, both of date) Can be used 
together with VHP and DAP to compute an 
initial flyby trajectory state, bu ? requires /(-plane 
aim point informaFon, e.g., B and 6 (deg). 

ZAPS Angle between arrival Y m vector and the arrival 
planet-to-Sun vector. Fquivalent to planet-probe- 
Sun angle at far encounter; for subsolar impact 
would be equal to 1K0 deg. Can he used with 
FTSP, VHP. DAP, and 6 to determine solar phase 
angle al periapse, entry, etc. (deg) 

ZAPF Angle between arrival I,. vector and the planet-to- 
Farth vector. Fquivalent to planet-probe-Ear'h 
angie at far encounter (deg). 

FTSP Angle in arrival /(-plane, measured tmm T-axis*, 
clockwise to projection of Sun-to-planct vector. 
Equivalent to solar occultation region centerline 
direction in /(-plane (deg). 

FTFP Angle in arrival W-plane. measured from T-axis , 
clockwise, to projection of harih-to-planet vector. 
Equivalent to Farth occultation region centerline 
direction in /( plane (deg). 

*HS? 4r.d LThP plot* are bated on T a a*, defined * be inn parallel 10 
elliptic plane (u.* text tor cxpUruttorii 
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V. Table of Constants 

('oi.»I;iiiIn used in pcth'iate ihc minimal ion ptesenled ate 
suiuuian/rd itt iIns section. 


A. Sun 

V#H 1,7 km 


knv 1 s ; 


B. Larth/Moon Systam 



-to.i Ml 5 >1 km’ s-' 

t '\i 

/ IK/7/ 

5»X <>tX) 44K O ’« km* s’ 

* 

0 0010S'<,.» 

f* 

\VI7ff 1 k'I 

i..rs 140 km 

Jupiter System 

r 'if 

'■ M.V/7M 

' 1 > ' 1; <*IK M4 kn» V 


’ 0.014 7 .UOOS 

" .11 rrn k 



«.\V. 

1 »1«*^ ■** 

• » 

Suit.H'C 
Kail his. 
km 

Mean 

Oiluf 

RjdmsJ** 

1 «1 
hua 

Period. 
d*i> .> 

Jupitei 

i :a fiHf* 46: 

7 1 .WH 0 


0.41 Win ' { 

Kilims 

$0 54 

is:o 

4 21 f* 71 

\ vh.ik: 

1 lliop.1 

514(* 

l.VU 

670<*88 

5 $M 181^ 

tfAinmeile 

*»88s 

>os 

1070J38 

7 1*4*5.10 

tjllMo 

7 i v 

:44S 

188:57*. 

16 68*W8 


(•H ompiilcil I tom period ami Jupiter t»W, minuted 
(••Hoi System III Jupitet lotation ule 

Direction of the Jovian planetary equatorial north pole 
(in l ; arth Mean Fquatoi of 1^50.0 coordinates): 

o # , ^.OOl^.ilcg 
~ M.5040 1 *. dejs 

0. Sources 

The constants represent the D1--1 IS planetary ephemeris 
(l^Sl) and the Voyager-2 Jupiter encounter trajectory 
reconstruction data ( 1070 ). 
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